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Abstract. We evaluate the prospects for detecting a non-standard light Higgs boson with a significant branch- 
ing ratio to two photons, in the Run II of the Fermilab Tevatron. We derive the reach for several channels: 27 
\ inclusive, 2y + 1 jet, and 27 -I- 2 jets. We present the expected Run II limits on the branching ratio of ft ^ 77 

, as a function of the Higgs mass, for the case of "bosonic", as well as "topcolor" Higgs bosons. 
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^ ■ MOTIVATION 

J> , The Standard Model (SM) is very economical in the sense that the Higgs doublet responsible for electroweak 

' symmetry breaking can also be used to generate fermion masses. The Higgs boson couplings to the gauge bosons, 
^5 I quarks, and leptons are therefore predicted in the Standard Model, where one expects the Higgs boson to decay 
mostly to b-jets and tau pairs (for low Higgs masses, Mh ^ 140 GeV), or to WW or ZZ pairs, (for higher Higgs 
I masses, Alh > 140 GeV). Since the Higgs boson is neutral and does not couple to photons at tree level, the branching 
■ ratio B(/i —> 77) is predicted to be very small in the SM, on the order of 10^"^ — 10^"*. 

I In a more general framework, however, where different sectors of the theory are responsible for the physics of 

. flavor and electroweak symmetry breaking, one may expect deviations from the SM predictions, which may lead to 
' drastic changes in the Higgs boson discovery signatures. One such example is the so called "fermiophobic" (also 
known as "bosophilic" or "bosonic") Higgs, which has suppressed couplings to all fermions. It may arise in a variety 
of models, see e.g. [1]. A variation on this theme is the Higgs in certain topcolor models, which may couple to heavy 
(■H I quarks only [2] . Some even more exotic possibilities have been suggested in the context of theories with large extra 
dimensions [3]. Finally, in the minimal supersymmetric standard model (MSSM), the width into bb pairs can be 
suppressed due to 1-loop SUSY corrections, thus enhancing the branching ratios of a light Higgs into more exotic 
signatures [4,5]. In all these cases, the Higgs boson decays to photon pairs are mediated through a.W 01 heavy quark 



^ I loop and dominate for Mh < 100 GeV [6]. In the range 100 < ^ 160, they compete with the WW* mode, while 
for Mh ^ 160 GeV, h WW completely takes over. Current bounds from LEP [7] are limited by the kinematic 
reach of the machine. The existing Run I analyses at the Tevatron have utilized the diphoton plus 2 jets [8-10] and 
inclusive diphoton [10] channels and were limited by statistics. Since they only looked for a "bosonic" Higgs [1], they 
did not consider the Higgs production mechanism through gluon fusion, which can be a major additional source of 
signal in certain models [2]. Since /i ^ 77 is a very clean signature, it will allow the Tevatron to extend significantly 
those limits in its next run. 

In this study we shall evaluate the Higgs discovery potential of the upcoming Tevatron runs for several diphoton 
channels. We shall concentrate on the following two questions. First, what is the absolute reach in Higgs mass as 
a function of the h — > 77 branching ratio? Second, which signature (inclusive diphotons, diphotons plus one jet, or 
diphotons plus two jets) provides the best reach. We believe that neither of those two questions has been adequately 
addressed in the literature previously. 



TEVATRON REACH FOR A BOSONIC HIGGS 



Here we consider the case of a "bosonic" Higgs, i.e. models where the Higgs couplings to all fermions are suppressed. 
Then, the main Higgs production modes at the Tevatron are associated Wh/Zh production, as well as WW/ZZ 
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FIGURE 1. The total background in the inclusive diphoton channel, as well as the individual contributions from 77, 7j 
and jj production. 



fusion. All of these processes have comparable rates [11], so it makes sense to consider an inclusive signature first 



Inclusive channel: analysis cuts 

We use the following cuts for our inclusive study: two photons with pxil) > 20 GeV and rapidity |r;(7)| < 2, 
motivated by the acceptance of the CDF or D0 detectors in Run II. Triggering on such a signature is trivial; both 
collaborations will have diphoton triggers that are nearly fully efficient with such offline cuts. 

We assume 80% diphoton identification efflciency, which we apply to both the signal and background estimates 
on top of the kinematic and geometrical acceptance. Again, this efficiency is motivated by the CDF/D0 EM ID 
efficiency in Run I and is not likely to change in Run II. 

Inclusive channel: background 

The main backgrounds to the inclusive diphoton channel come from the QCD production of dijets, direct photons, 
and diphotons. In the former two cases a jet mimics a photon by fragmenting into a leading tt^ /r] meson that further 
decays into a pair of photons, not resolved in the calorimeter. 

We used the PYTHIA [12] event generator and the experimentally measured probability of a jet to fake a photon [8] 
to calculate all three components of the QCD background. The faking probability depends significantly on the 
particular photon ID cuts, especially on the photon isolation requirement (see, e.g. [8,13,14]). For this study we used 
an E'x-dependent jet-faking-photon probability of 



which is obtained by taking the ?7-averaged faking probabilities used in the D0 Run I searches [8]. The fractional 
error on P(jet 7) is about 25% and is dominated by the uncertainty on the direct photon fraction in the jet + 7 
sample used for its determination. (For high photon Et, however, the error is dominated by the available statistics.) 



[10]. 





This probability is expected to remain approximately the same in Run II for both the CDF and D0 detectors. We 
used 80% ID efficiency for the pair of photons, and required the photons to be isolated from possible extra jets in 
the event. We accounted for NLO corrections via a constant fc-factor of 1.34. 

Adding all background contributions, for the total background in the inclusive diphoton channel we obtain the 
following parametrization: 

+ (rSv) + ( ] '""^ + ( } ' 

where pi = 6.45, ps = -0.029, = 2.44, = -0.011 and ps = 0.00005. In the region M^^ > 100 GeV it is 
dominated by direct diphoton production and hence is irreducible. The expected statistical plus systematic error on 
this background determination is at the level of 25%, based on the jet-faking photon probability uncertainty. For 
larger invariant masses, however, the accuracy is dominated by the uncertainties in the direct diphoton production 
cross section, which will be difficult to measure independently in Run II, so one will still have to rely on the NLO 
predictions. On the other hand, for narrow resonance searches one could do self-calibration of the background by 
calculating the expected background under the signal peak via interpolation of the measured diphoton mass spectrum 
between the regions just below and just above the assumed resonance mass. Therefore, in our case the background 
error will be purely dominated by the background statistics. A combination of the interpolation tecnique and the 
shape information from the theoretical NLO calculations of the direct diphoton cross section is expected to result 
in significantly smaller background error in Run II. 

The total background, as well as the individual contributions from 77, jj and jj production, are shown in 
Fig. 1. Additional SM background sources to the inclusive diphoton channel include Drell-Yan production with both 
electrons misidentified as photons, W77 production, etc. and are all negligible compared to the QCD background. 
The absolute normalization of the background obtained by the above method agrees well with the actual background 
measured by CDF and D0 in the diphoton mode [10,14]. 

In Fig. 2 we show the 95% CL upper limit on the differential cross section after cuts d{e x (7(77 + X))/dM^^ as 
a function of the diphoton invariant mass M^^, given the above background prediction (here e is the product of the 
acceptance and all efficiencies). This limit represents 1.96cr sensitivity to a narrow signal when doing a counting 
experiment in 1 GeV diphoton mass bins. This plot can be used to obtain the sensitivity to any resonance decaying 
into two photons as follows. One first fixes the width of the mass window around the signal peak which is used in 
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FIGURE 3. Significance S/\fB (in arbitrary units), as a function of the mass window width w (in units of F), for 
Breit-Wigner or Gaussian resonance. 



the analysis. Then one takes the average value of the 95% C.L. limit in da/dMjj across the mass window from 
Fig. 2 and multiplies it by ^/w/GeV, where w is the width of the mass window^ , to obtain the corresponding 95% 
CL upper limit on the signal cross-section after cuts. Similar scaling could be used if one is interested in the 3(T or 
5cr reach. 



What is the optimum mass window cut? 



When searching for narrow resonances in the presence of large backgrounds (J5), the best sensitivity toward signal 
{S) is achieved by performing an unbinned maximum likelihood fit to the sum of the expected signal and background 
shapes. However, simple counting experiments give similar sensitivity if the size of the signal "window" is optimized. 
For narrow resonances the observed width^ F is dominated by the instrumental effects, and is often Gaussian. The 
background in a narrow window centered on the assumed position Mq of the peak in the signal invariant mass 
distribution could be treated as linear. Therefore, the Gaussian significance of the signal, S/^/B, as a function of 
the window width, w, is given by: 
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The function (1) is shown in Fig. 3, and has a maximum at w ~ 1.2F, which corresponds to ±1.2(r/2) cut around 
the resonance maximum. 

For resonances significantly wider than the experimental resoluton, the shape is given by the Breit-Wigner function, 
and in this case the significance is: 
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The square root enters the calculation since the significance is proportional to the background to the —1/2 power. 
Notice that the width is defined so that the cross-section at ±r/2 away from the peak is a factor of 2 smaller than the 
peak value (FWHM). For a Gaussian resonance the width is related to the variance a by F = 2 cr \/ In 4 ~ 2.35a. 



This function, also shown in Fig. 3, peaks at a similar value of w (w ~ 1.4r). We see that for both Gaussian and 

Breit-Wigner resonances, the significance does not appreciably change when using a w = IF — 2r cuts. For our 
analysis we shall use two representative choices: w — 1.2F and w — 2T for the mass window, which we shall always 
center on the actual Higgs mass. 

Clearly, one can do even better in principle, by suitably resizing and repositioning the mass window around the 
bump in the combined S + B distribution. Because of the steeply falling parton luminosities, the signal mass peak 
is skewed and its maximiun will appear somewhat below the actual physical mass. In our analysis we choose not to 
take advantage of these slight improvements, thus accounting for unknown systematics. 

Inclusive channel: results 

In Tables 1 and 2 wc show the inclusive jj + X background rates in fb for different Higgs masses, for w = 1.2F and 
w = 2r mass window cuts, respectively. Here we have added the intrinsic width F^ and the experimental resolution 



TABLE 1. Background rates in fb for w = 1.2r mass cut, and significance {S/^/B, for 1 fb ^ of data, and assuming 
B{h — > 77) = 100%) as a function of the Higgs mass M^. The signal consists of associated Wh/Zh production and 
WW/ZZ fusion. 



(GeV) 


77 + X 
bknd 
(fb) 


Significance S/Vb 


■y-f + X 


77+1 jet 


77 + 2 jots 


PT > 20 


PT > 25 


PT > 30 


Pt > 35 


PT > 20 


PT > 25 


Pt > 30 


PT > 35 


100. 


271.7 


16.5 


31.3 


34.2 


36.3 


35.4 


31.9 


36.4 


35.1 


31.1 


120. 


166.6 


13.0 


24.4 


26.7 


28.5 


28.5 


24.7 


29.5 


29.2 


26.7 


140. 


103.0 


10.7 


20.1 


22.4 


23.9 


23.9 


20.6 


24.0 


23.7 


21.0 


160. 


64.3 


8.9 


17.0 


19.1 


20.1 


20.4 


17.2 


20.2 


20.8 


19.2 


180. 


40.5 


7.3 


13.5 


15.0 


16.2 


16.5 


13.6 


16.3 


16.7 


16.2 


200. 


26.1 


5.8 


10.6 


11.9 


12.5 


12.7 


10.4 


12.3 


12.7 


12.0 


250. 


9.4 


3.7 


6.7 


7.5 


7.9 


8.2 


6.6 


7.9 


8.7 


8.4 


300. 


4.8 


2.2 


3.8 


4.3 


4.7 


4.7 


3.6 


4.2 


4.9 


4.7 


350. 


2.3 


1.5 


2.8 


3.1 


3.3 


3.3 


2.4 


2.8 


3.0 


3.4 


400. 


i.2 


1.0 


i.cS 


2.0 


2.0 


2.1 


1.7 


1.7 


2.1 


2.1 



Fexp = 2Vhi4 X (Tcxp ^ 2.35 x 0.15V2^/E{j) ~ in quadrature: F = (F^ + F2^p)^''^ The width F varies 

between 3.5 GeV for Mh = 100 GcV and 29.0 GeV for = 400 GeV. The two tables also show the significance (for 
1 fb~^ of data, and assuming B{h 77) = 100%) in the inclusive diphoton channel when only associated Wh/Zh 
production and WW/ZZ h fusion are included in the signal sample. We see that (as can also be anticipated from 
Fig. 3) a, w = 1.2F cut around the Higgs mass typically gives a better statistical significance, especially for lighter 
(and therefore more narrow) Higgs bosons. 

Exclusive channels: analysis 

The next question is whether the sensitivity can be further improved by requiring additional objects in the event. 
The point is that a significant fraction of the signal events from both associated Wh/Zh production and WW/ZZ 
fusion will have additional hard objects, most often QCD jets. In Fig. 4 we show the "jet" multiplicity in associated 
Wh production, where for detector simulation we have used the SHW package [15] with a few modifications as in 
[16]. Here we treat "jets" in a broader context, including electrons and tau jets as well. 

Previous studies [10,9] have required two or more additional QCD jets. Here wc shall also consider the signature 
with at least one additional "jet", where a "jet" is an object with \ri\ < 2. The advantages of not requiring a second 
"jet" are twofold. First, in this way we can also pick up signal from WW/ZZ — > h fusion, whose cross-section 
does not fall off as steeply with Mh, and in fact for > 200 GeV is larger than the cross-section for associated 



TABLE 2. The same as Table 1, but for a w = 2T mass window. 



Mh 


77 + X 
bknd 

(fb) 


Significance S/^/B 


77 + X 


77 + 1 jet 


77 + 2 jots 


PI > 20 


PI > 25 


pi > ■'!() 


Pi > '■')!'> 


PI > 21) 


Pi > 2.") 


PT > :!() 


pi' > .35 


100. 


453.4 


UA 


27.0 


29.8 


31.7 


30.6 


27.7 


31.9 


30.5 


26.4 


120. 


278.1 


11.3 


21.3 


23.5 


25.1 


24.9 


21.9 


25.5 


25.1 


22.5 


140. 


171.9 


9.3 


17.5 


19.5 


21.0 


21.0 


17.7 


20.7 


21.0 


18.5 


160. 


107.3 


8.0 


15.1 


16.7 


18.0 


18.2 


15.4 


17.8 


18.2 


17.3 


180. 


67.6 


6.6 


12.2 


13.7 


14.6 


14.9 


12.4 


14.3 


15.2 


14.2 


200. 


43.6 


5.4 


10.1 


11.4 


12.1 


12.1 


9.9 


11.9 


12.3 


11.1 


250. 


15.7 


3.6 


6.5 


7.3 


7.7 


7.8 


6.4 


7.6 


8.5 


8.0 


300. 


8.1 


2.1 


3.8 


4.2 


4.6 


4.5 


3.6 


4.3 


4.5 


4.5 


350. 


3.9 


1.6 


2.7 


3.0 


3.4 


3.4 


2.6 


3.1 


3.5 


3.1 


400. 


2.1 


1.1 


1.9 


2.2 


2.3 


2.4 


1.6 


2.1 


2.6 


2.4 



Wh/Zh production"^. Events from WW/ZZ h fusion typically contain two very hard forward jets, one of which 
may easily pass the jet selection cuts. In Fig. 5 we show the pseudorapidity distribution of the two spectator jets in 
WW/ZZ — > h fusion (red) and associated Wh/Zh production (blue). Second, by requiring only one additional jet, 
we win in signal acceptance. In order to compensate for the corresponding background increase, we shall consider 
several px thresholds for the additional jet, and choose the one giving the largest significance. 

For the exclusive channels we need to rescale the background from Fig. 1 as follows. From Monte Carlo we obtain 
reduction factors of 4.6 ± 0.5, 6.2 ± 1.0, 7.6 ± 1.4, and 8.6 ± 1.5 for the 77 + 1 jet channel, with prij) > 20, 25, 30 
and 35 GeV, respectively. For the 77+2 jets channel the corresponding background reduction is 21 ±5, 38 ± 12, 
58 ± 21, and 74 ± 26, depending on the jet pt cuts. These scaling factors agree well with those from the CDF and 
D0 data from Run I. 

Notice that we choose not to impose an invariant dijet mass {Mjj) cut for the 77 + 2 jets channel. We do not 
expect that it would lead to a gain in significance for several reasons. First, given the relatively high jet pT cuts 
needed for the background suppression, there will be hardly any background events left with dijet invariant masses 
below the (very wide) W/Z mass window. Second, the signal events from WW/ZZ fusion, which typically comprise 
about 25 — 30% of our signal, will have a dijet invariant mass distribution very similar to that of the background. 
Finally, not imposing the Mjj cut allows for a higher signal acceptance because of the inevitable combinatorial 
ambiguity for the events with > 2 jets. 

The significances for the two exclusive channels, with the four different jet px cuts, are also shown in Tables 1 
and 2. We see that the exclusive 77 + 2 jets channel with pxij) > 30 GeV typically gives the largest significance, 
but our new exclusive 77 + 1 jet channel is following very close behind. 

Exclusive channels: results 

We are now ready to present our results for the Run II Tevatron reach for a bosonic Higgs. In Fig. 6 we show 
the 95% CL upper limit on the branching ratio B{h 77), with 0.1 (cyan), 2.0 (green) and 30 fb~^ (red), as a 
function of M/j. For each mass point, we compare the significance for both the inclusive as well as the exclusive 
channels with all the different cuts, and for the limit we choose the channel with the set of cuts providing the best 
reach. It turns out that for the case at hand the winners are: o: 27 + 2j, with priJ) > 25 GeV; □: 27 + 2j, with 
PtU) > 30 GeV, and O: 27 + Ij, with prij) > 30 GeV. In the figure we also show the HDECAY [17] prediction 
for B(/i — > 77) in case of a "bosonic" Higgs. The reach shown for 0.1 fb~^ is intended as a comparison to Run I, 
in fact for the 0.1 fb~^ curve we scaled down both the signal and background cross-sections to their values at 1.8 
TeV center-of-mass energy, keeping the efficiencies the same. In other words, the region marked as Run F would 
have been the hypothetical reach in Run I, if the improved Run II detectors were available at that time. As seen 



' In the case of a topcolor Higgs (see the next section) we would also pick up events with initial state gluon radiation, 
comprising about 30% of the gluon fusion signal, which is the dominant production process for any Higgs mass. 
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FIGURE 4. The number of "jets" , which stands for QCD jets, tan jets and electrons, in associated Wh production, once 
we require the two photons from the Higgs to pass the photon ID cuts. 



from Fig. 6, the reach for a "bosonic" Higgs [1] (at 95% CL) in Run Ila and Run lib is ~ 115 GeV and ~ 125 GeV, 
correspondingly. This is a significant improvement over the ultimate reach from LEP [7] of 105 GeV. 

TEVATRON REACH FOR A TOPCOLOR HIGGS 

Here we consider the case of a "topcolor" bosonic Higgs, where the Higgs also couples to the top and other heavy 
quarks [2] . We therefore include events from gluon fusion into our signal sample. We used the next-to-leading order 
cross-sections for gluon fusion from the HIGLU program [18]. 

In Tables 3 and 4 wo show the significance (for 1 fb~^ of data, and again assuming B(/i 77) — 100%) in the 
inclusive and the two exclusive channels, for the topcolor Higgs case. Since gluon fusion, which rarely has additional 
hard jets, is the dominant production process, the inclusive channel typically provides the best reach. However, 
the 27 + Ij channel is again very competitive, since the additional hard jet requirement manages to suppress the 
background at a reasonable signal cost. We see that our new 27 + Ij channel clearly gives a better reach than the 
27 + 2j channel [8-10]. For Higgs masses above ~ 180 GeV, it sometimes becomes marginally better even than 
the inclusive diphoton channel. The specific jet pt cut and mass window size w seem to be less of an issue - from 
Tables 3 and 4 we see that prij) > 25, prij) > 30 GeV and pxij) > 35 GeV work almost equally well, and for 
Mfi > 200 GeV both values of w are acceptable. 

In Fig. 7 we show the Run II reach for the branching ratio B(ft, — !■ 77) as a function of the Higgs mass, for the 
case of a "topcolor" Higgs boson. This time the channels with the best signal-to-noise ratio are: o: inclusive 27 + X, 
and □: 27 -|- Ij, with priJ) > 30 GeV; both with w = 1.2r. 

CONCLUSIONS 

We have studied the Tevatron reach for Higgs bosons decaying into photon pairs. For purely "bosonic" Higgses, 
which only couple to gauge bosons, the 27 + 2j channel offers the best reach, but the 27 + Ij channel is almost as 
good. For topcolor Higgs bosons, which can also be produced via gluon fusion, the inclusive 27 -|- X channel is the 
best, but the 27 -|- Ij channel is again very competitive. We see that in both cases the 27 -|- Ij channel is a no-lose 
option! 
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FIGURE 5. Pseudorapidity distribution of the two spectator jets in WW/ZZ h fusion (red) and associated Wh/Zh 
production (blue). The boxed region represents the off-hne selection cuts. 
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TABLE 3. The same as Table 1, but for a topcolor Higgs, i.e. gluon fusion events are included in the signal. 



Mh 
(GeV) 


77 + X 
bknd 
(fb) 


Significance S/\/B 


77 + X 


- - jot 


— 2 ji'ts 


PT > 20 


PT > 25 


PT > 30 


PT > 35 


PT > 20 


PT > 25 


PT > 30 


PT > 35 


100. 


271.7 


54.3 


43.4 


44.6 


45.1 


42.5 


33.6 


37.8 


36.4 


32.1 


120. 


166.6 


42.5 


35.3 


36.9 


37.3 


35.3 


26.5 


30.8 


30.2 


27.5 


140. 


103.0 


35.5 


30.4 


32.1 


33.0 


31.6 


22.3 


25.6 


25.1 


22.0 


160. 


64.3 


30.1 


27.3 


28.8 


29.2 


28.3 


19.5 


22.4 


22.4 


20.3 


180. 


40.5 


24.5 


22.4 


23.8 


24.5 


23.8 


15.7 


18.2 


18.3 


17.6 


200. 


26.1 


19.9 


18.6 


19.9 


20.3 


20.0 


12.7 


UA 


14.3 


13.3 


250. 


9.4 


12.9 


12.5 


13.3 


13.6 


13.4 


8.4 


9.7 


10.2 


9.8 


300. 


4.8 


9.1 


8.7 


9.3 


9.6 


9.4 


4.9 


5.6 


5.9 


5.9 


350. 


2.3 


8.4 


7.9 


8.4 


8.4 


7.9 


4.0 


4.1 


4.5 


4.5 


400. 


1.2 


6.4 


6.1 


6.1 


6.3 


6.2 


3.0 


3.0 


3.5 


3.1 



TABLE 4. The same as Table 3, but for a w = 2r mass window. 



Mh 
(GeV) 


77 + X 
bknd 
(fb) 


SigiiificaiKx' S/s^ti 


77 + X 


77 + 1 jet 


77 + 2 jets 


PT > 20 


PT > 25 


PT > 30 


PT > 35 


PT > 20 


PT > 25 


PT > 30 


PT > 35 


100. 


453.4 


47.7 


37.6 


38.9 


39.5 


36.9 


29.0 


33.0 


31.6 


27.1 


120. 


278.1 


37.7 


31.2 


32.6 


33.1 


31.2 


23.5 


26.6 


25.9 


23.1 


140. 


171.9 


31.4 


26.7 


28.2 


29.0 


27.9 


19.2 


22.1 


22.0 


19.4 


160. 


107.3 


26.6 


24.3 


25.3 


26.0 


25.2 


17.4 


19.7 


19.5 


18.3 


180. 


67.6 


22.3 


20.5 


21.7 


22.1 


21.6 


14.2 


16.2 


16.8 


15.4 


200. 


43.6 


18.8 


17.7 


19.1 


19.4 


18.9 


12.0 


13.7 


13.8 


12.5 


250. 


15.7 


12.8 


12.1 


13.2 


13.4 


13.1 


8.2 


9.3 


10.0 


9.6 


300. 


8.1 


9.0 


8.5 


9.1 


9.2 


9.0 


5.0 


5.6 


5.6 


5.7 


350. 


3.9 


8.3 


7.6 


8.1 


8.1 


7.9 


4.0 


4.7 


4.6 


4.4 


400. 


2.1 


6.3 


6.1 


6.3 


6.3 


6.1 


2.9 


3.4 


3.7 


3.6 




FIGURE 7. The same as Fig. 6, but for a topcolor Higgs, i.e. gluon fusion events are included in the signal. The channels 
with the best S/y/B ratio are: o: inclusive 27 + X, and □: 27 + Ij, with prU) > 30 GeV; both with w = 1.2 F. 



